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ABSTRACT 



We compare the predictions of simple models for the radio emission from classical 
*-£h ■ novae with the MERLIN radio observations of nova V723 Cas (Nova Cas 1995). Spher- 

ically symmetric and ellipsoidal radiative transfer models are implemented in order to 
Q . generate synthetic emission maps. These are then convolved with an accurate repre- 

ss ' sentation of the uv coverage of MERLIN. The parameters and geometry of the shell 

c/3 \ model are based on those returned by fitting models to the observed light curve. This 

allows direct comparison of the model images with the nine 6 cm MERLIN images of 
V723 Cas. 

It is found that the seemingly complex structure (clumping, apparent rotation) 
■ evident in the observations can actually be reproduced with a simple spherical emission 

[■ — ' model. The simulations show that a 24 hour track greatly reduces the instrumental 

r**"* . effects and the synthetic radio map is a closer representation of the true (model) sky 

brightness distribution. It is clear that interferometric arrays with sparse uv coverage 
(e.g. MERLIN, VLBA) will be more prone to these instrumental effects especially 
OO ■ when imaging ring-like objects with time dependent structure variations. A modelling 

' approach such as that adopted here is essential when interpreting observations. 

o. ; 
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1 INTRODUCTION et al, 1996, 2000, 2005; Heywood et al, 2002, 2005). Both 

radio and optical images of classical novae often show non- 
Classical novae are the result of thermonuclear explosions in homogeneous, albeit spherical or ellipsoidal structures. A 
interacting binary stars. The system consists of a hot white feature common to all MERLIN observations of novae is ap- 
dwarf primary with a main sequence secondary companion. parent rotation between epochs. 

The secondary is thought to fill its Roche lobe and trans- In thig paper we investigate the struc ture seen in MER- 
fers hydrogen-rich material via an accretion disc onto the LIN observations of V723 Cas (Nova Cas 1995) using models 
surface of the white dwarf. When the pressure in the de- fitted tQ the radio light curve (Heywood et al 20 05 - here- 
generate material at the base of the accreted layer reaches after p aper T) The simulations present ed in this paper pro- 
a critical value a thermonuclear runaway takes place (Star- vide an ingight into how the radio imageg that regult from 
rfield 1989), ejecting -10 -10 M of material into the interferometer observations relate to the true distribution of 
interstellar medium at velocities of up to several thousand the emitting materia i. Radio images are often cited during 
kilometres per second. This causes an increase m the bolo- discussions of the morphology of nova shells. Although it is 
metric luminosity of the system of a few orders of magnitude common knowledge that interferometer images are an ap- 
over a timescale of a few days. proximation of the true sky brightness distribution, it has 
Radio emission from the nova starts to develop at a later not previously been determined just how many of the fea- 
time than the optical. The white dwarf star photoionises tures which are common to radio observations of novae are 
the ejecta giving rise to free-free emission at temperatures due to the interferometer and the aperture synthesis pro- 
of approximately 10 4 K. The evolution of the ejecta is gen- cess. The goal for this investigation is to generate synthetic 
erally rapid, allowing radio interferometric observations to emission maps based on shell models to be subsequently con- 
track the expanding shell as it brightens and then fades volved with an accurate uv coverage model of the MERLIN 
over favourable timescales (e.g. Taylor et al, 1987; Eyres interferometer. 
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Descriptions of the shell models used are presented in 
Section 2. A radiative transfer solution is implemented for 
both a spherically symmetric and an ellipsoidal shell based 
on the models described in Paper I. This process results 
in simulated emission maps which are then convolved with 
a uv model. In all cases this process, described in Section 
3, is consistent with the genuine MERLIN observations of 
V723 Cas presented in Paper I and reproduced in Figure 
1. This facilitates direct comparisons between the real and 
simulated radio images. 



2 SIMULATIONS 

Model images of the sky brightness distribution are gen- 
erated by solving the line of sight equations of radiative 
transfer through a model emitting region. Descriptions of 
the radiative transfer code and the model shell geometries 
are presented in this section, together with images showing 
the predicted brightness distributions. Synthetic interferom- 
eter maps are derived by masking the Fourier transform of 
these images with the uv coverage of the interferometer and 
then inverting the transform, this process is described in 
Section 3. 




Figure 2. Coordinate system for the axisymmetric radiative 
transfer code. This figure shows how the shell is defined in the 
cylindrical polar coordinates used by the model. 



2.1 Axisymmetric radiative transfer 

The radiative transfer code uses cylindrical polar coordi- 
nates to define the emitting region as used by, for example, 
Dougherty et al (2003). The coordinate system is shown in 
Figure 2. The shell model is constructed by creating a den- 
sity function in the plane defined by R > and z > 0. A 
3D shell is created by mirroring the distribution about the 
plane defined by R and <f> and rotating through 360 degrees 
about the z axis. Thus the quarter circles indicated by the 
solid lines in Figure 2 would define full spherical regions. 
The configuration shown in Figure 2 is representative of the 
Hubble flow model described in Paper I (and previously by 
Seaquist & Palimaka, 1977, and by Hjellming et al, 1979) 
and is discussed further below. 

The density and temperature are defined in an array of 
dimensions 256 x 512 in R and z respectively and when the 
coordinate transform is applied to define the 3D emitting re- 
gion the result is mapped on to a 256 x 256 pixel image. The 
arrays defining the image are chosen to be 2™ in size, where 
n is an integer, so that the Fast Fourier Transform (FFT) 
technique can be applied when convolving these images in 
the uv plane. 

The frequency of the radiation is specified and the radia- 
tive transfer equations are solved along gridded lines of sight 
through the shell. The code then calls a Hierarchical Data 
System (HDS) library which outputs images of intensity, op- 
tical depth and brightness temperature in the STARLINK 
N-dimensional Data Format (NDF). 

An inclination angle for the orientation of the model 
can also be specified. This is denoted by 6 and is the angle 
between the v axis in the plane of the sky and the z axis of 
the coordinate system. 



2.2 Spherically symmetric Hubble flow model 

The Hubble flow model applied to V723 Cas in Paper I is 
used as the basis for the first run of the axisymmetric radia- 
tive transfer code. The values of mass and inner and outer 
ejection velocity returned by the best fit routine are used to 
define the density and the shell geometry for this model. The 
inner and outer shell boundaries, and consequently the shell 
density, are time dependent and emission maps are created 
with epochs corresponding to the 6 cm MERLIN observa- 
tions presented in Paper I. 

The density function of the Hubble flow model is given 

by 

p{r) = 4^^T) (1) 

which is a function of r, the radial distance from the centre 
of the shell, and also of n and r 2 , the radii of the inner and 
outer shell boundaries respectively. Hence when the shell is 
defined in the radiative transfer code the density array is 
filled according to 



p(R,z) 



M 



(2) 



4tt(R 2 + z 2 )(r 2 -ri) 

for n < r < T2- This is a function of R and z and the 
parameters r 2 and n which are fixed for a given epoch. 
Similarly for n < r < r 2 the temperature is fixed at 17000 K 
(see Paper I). Outside the shell boundaries the density is 
assumed to be zero and the temperature is set at 100 K. The 
inclination angle is a redundant parameter in the spherical 
case. 

The nine model intensity maps created by this process 
are presented in Figure 3 where it can be seen how initially 
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Figure 1. 6 cm MERLIN images of the nova V723 Cas, as published by Heywood ct al. (2005), reproduced here for comparison with 
the simulations. The contour levels are (-3, 3, 3^2, 6, 6\/2, 12, 12\/2...) X 0.15 mjy. 
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the emission is entirely optically thick. As the material ex- 
pands the outer layers start to become optically thin and 
eventually the inner boundary of the shell is revealed as a 
bright ring. 



2.3 Ellipsoidal shell model 

Some observations of classical novae show prolate ellipsoidal 
shells (e.g. FH Ser, Gill & O'Brien, 2000) and the intention 
here is to determine whether it is possible to distinguish 
between a spherical shell or an ellipsoidal shell by studying 
the radio images. 

The shell model is based on the spherical Hubble flow 
model as described in Section 2.2 however some additional 
parameters are introduced. A polar and equatorial direction 
are defined which, with reference to Figure 2, are along the 
z and R axes respectively. The polar velocity, v p , is greater 
than the equatorial velocity, v e , and the two are related by 
the factor V such that v e = V v p . The velocity in these 
regions is a function of the angle ip which is defined as 



ip = tan 



R 



at ip = and v(ip) = v p at ip 



(3) 
tt/2, 



such that v(tp) ■ 
see Figure 4. 

The velocity as a function of ip is governed by the ellipse 
equations. In the general case, for an ellipse of semi-major 
axis a, semi-minor axis b and eccentricity e, the radius r 
as a function of 9, the angle between the semi-minor and 
semi-major axes, is given by 



aVl - e 2 
VI -e 2 cos 2 9 
where 



(4) 



(5) 



By subtituting Equation 5 into Equation 4 it can be shown 
that 



b 2 



ab 



V^ 2 



b 2 ) cos 2 



(6) 



and the above formula is used to govern how the velocity 
varies as a function of ip in the shell, i.e. 



v(ip) 



\Jvl — (yl — v 2 ,) cos 2 ip 



(7) 



The elliptical shell models are again constructed based 
on the best fit parameters of V723 Cas as returned by the 
light curve fitting in Paper I. The best fitting inner and outer 
shell velocities are taken to be the velocities of the inner and 
outer shell boundaries at the pole. The shell mass is taken to 
be the best fitting shell mass from the spherical model. The 
value of V is arbitrarily chosen to be 0.6 and the inclination 
of the shell to the plane of the sky is taken to be zero degrees 
i.e. the shell is viewed 'side-on'. 

As with the spherically symmetric model, nine simu- 
lated maps are calculated to correspond to the nine epochs 
of 6 cm MERLIN observations. Equation 7 is used to deter- 
mine vi(ip) and V2(tp), the velocities of the inner and outer 
shell boundaries as a function of ip, and these velocities are 




Figure 4. Geometry for elliptical shell model. 

then used to calculate n(ip) and r2(ip), the radii of the in- 
ner and outer shell boundaries for the specific epoch. The 
density function 

M 



p{R,z) 



(8) 



4tt(R 2 + z 2 )(r 2 -n) 

is then applied to populate the density array in R and z. 
Once again the temperature is set at 17000 K within the 
shell boundaries and is assumed to be 100 K elsewhere. 

The nine model intensity distribution maps are pre- 
sented in Figure 5. As with the spherical case the emission is 
initially optically thick before expansion causes the optical 
depth to decrease and the inner shell boundary is revealed. 
An interesting feature to note is the 'waist' of brighter emis- 
sion in the equatorial region. This is due to the inverse de- 
pendence on the quantity r 2 — ri in the density equation. 
This is smaller for small values of ip than it is for values of 
ip close to 7r/2 radians giving rise to a density enhancement 
near the equator of the shell. 



3 SIMULATED MERLIN IMAGES 

Once the model images have been calculated the im- 
age files have to be correctly formatted before they can 
be imported into the National Radio Astronomy Obser- 
vatory's Astronomical Image Processing System (AIPS; 
http://www.aoc.nrao.edu/aips/) which is used to gener- 
ate the synthetic MERLIN images. A World Coordinate Sys- 
tem, frequency and Stokes parameter axes are added. The 
AIPS task UVCON* reads the properties of the array from 

* Recently a bug was discovered in the AIPS task UVCON 
whereby in certain cases visibilities in the simulated data set 
would be recorded with telescopes at negative elevations, and 
then conversely visibilities would be missing from high elevations 
where they should be present (L. Kogan, private communication). 

When this problem came to light the simulations in this pa- 
per were naturally re-examined to check for this problem. This 
footnote is merely to confirm that the results presented are not 
affected. The declination of V723 Cas coupled with the latitudes 
of the stations in the MERLIN array result in positive elevations 
throughout, even for a 24 hour track. 
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Figure 3. Model intensity distributions for the Hubble flow model. These maps are calculated at times following outburst shown in each 
panel, coinciding with the 6 cm MERLIN observations of V723 Cas. The greyscale is in units of mjy. 
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Figure 5. Model intensity distributions for the elliptical shell model, see Fig. 3. 
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an antenna file and determines the uv coverage of the inter- 
ferometer. The visibilities are then calculated from a gridded 
FFT of a model image. Since the models were defined by fits 
to the light curve, the flux density is very close to that mea- 
sured in the actual observations. Gaussian noise is added to 
the data, the level of which is calculated by UVCON based 
on integration times and the total length of the track. The 
uv data produced by UVCON are deconvolved, cleaned and 
imaged using IMAGR. In all cases where possible, the imag- 
ing process has been consistent with the methods used to 
produce the 6 cm MERLIN maps shown in Figure 1. The 
cleaning process was terminated at the first negative clean 
component and a circular restoring beam of 50 mas was 
used. 

Figure 6 shows the simulated MERLIN images based on 
the spherically symmetric Hubble flow model. These images 
are based on a 12 hour track which is of similar length to 
the original MERLIN observations. The images in Figure 7 
take the same model intensity distributions as those used to 
generate Figure 6 except the simulated uv coverage is based 
on a 24 hour track instead of a 12 hour track. Figure 8 shows 
the model MERLIN observations of the elliptical shell model 
based on a 12 hour track, and finally Figure 9 shows the 
same model with a 24 hour track. In each 3x3 grid, left to 
right, top to bottom, the images correspond to the December 
1996, January 1997, March 1998, December 1998, February 
2000, April 2000, January 2001, June 2001 and October 2001 
MERLIN observations of V723 Cas. Contour levels are (-3, 
3, 3\/2, 6, 6\/2, 12, 12y/2...) x 0.15 mjy for each image, as 
quoted in the captions. 



4 DISCUSSION 

Comparison of Figure 6 to the MERLIN images presented 
in Figure 1 reveals that even for this simple spherical shell 
model the similarities between the simulated radio observa- 
tions and the real MERLIN observations are striking. Ini- 
tially the source is point-like, not much bigger than the 
restoring beam. The source expands and as the turn over 
into optically thin emission begins apparent structure in the 
shell becomes evident. The source becomes ring-like then di- 
vides into the double-peaked structure which is common to 
all multi-epoch MERLIN observations of classical novae. 

The most intriguing feature of these simulated maps 
is the switching of the alignment of the dominant peaks in 
the map, giving the impression of rotation. In Paper I we 
speculated whether this was due to optical depth effects, 
with higher density regions of the shell remaining optically 
thick for longer and appearing as peaks of radio emission 
as the photosphere recedes. Here this explanation can be 
ruled out as the mass distribution supplied by the model is 
spherical with a smooth radial density profile. 

Comparing Figure 7 to Figure 6, shows that for a 24 
hour observation the radio emission in the synthetic maps 
is more circularly symmetric and consquently more faithful 
to the model images. The effect whereby the alignment of 
the dominant peaks in emission appears to rotate between 
epochs is also almost completely eradicated. The increase 
in uv coverage provided by the 24 hour observation must 
be the reason for this, providing both an improved signal 
to noise ratio and a drastically altered sidelobe structure in 



Table 1. Comparison of the background rms values in mjy/beam 
for both the 12 and 24 hour simulations of the spherical Hubble 
flow model. 



Epoch 


12 hour track 


24 hour track 


Dec 96 


0.169 


0.124 


Jan 97 


0.129 


0.119 


Mar 98 


0.177 


0.137 


Dec 98 


0.185 


0.141 


Feb 00 


0.141 


0.134 


Apr 00 


0.137 


0.133 


Jan 01 


0.136 


0.124 


Jun 01 


0.131 


0.117 


Oct 01 


0.128 


0.112 



the beam. The beam shape and uv coverage are discussed 
in Section 4.1. Table 1 lists the root mean square values of 
the background noise for both 12 and 24 hour tracks of the 
simulated spherical Hubble flow model. 

The simulated radio maps for the ellipsoidal shell model 
reveal a north-south bias with a double peaked source being 
present for the first five epochs. This gives way to a 'partial 
shell' structure which in turn is replaced by the east-west 
aligned peaks in the final map. Although a stronger north- 
south alignment of the emission would be expected from this 
ellipsoidal shell model the evolution of the radio emission 
over time displays many similiarities to the spherical case. 
Both initially develop a north-south bias which gives way 
to an east-west alignment and although this process occurs 
later for the ellipsoidal shell it is unlikely that a 12 hour 
track can be used to distinguish between an ellipsoidal or a 
spherical shell for objects of this nature. 

In order to determine whether a 24 hour track would 
produce a distinction between ellipsoidal and spherical shells 
the simulations presented in Figure 9 can be examined. For 
the optically thick stage the MERLIN images accurately rep- 
resent the simulated brightness distributions but for subse- 
quent epochs the images detect only the brighter east-west 
aligned emission around the inner shell boundary. The faint, 
extended structure along the major axis of the shell is be- 
low detectability. With a greater signal to noise ratio, or 
alternatively an object brighter than V723 Cas, it is likely 
that observations based on a 24 hour track could be used to 
examine the eccentricity of an ellipsoidal shell. 

4.1 Comparing the beams and the uv coverage 

Figure 10 shows the improved uv coverage of the 24 hour 
track compared to the 12 hour track. Plots of the beam 
shape for both 12 and 24 hour observations are presented 
in Figure 11. Immediately apparent is the effect that the in- 
creased uv coverage has on the sidelobes of the beam. The 
sidelobe structure in the 24 hour observation completely sur- 
rounds the central maximum and only deviates from circu- 
larity due to the declination of the source. The beam of 
the 12 hour simulation has S'-shaped sidelobe structure. As 
the source expands through the beam different parts of the 
shell will be buried in the sidelobes. The non-circular beam 
shape associated with the 12 hour track could be a possible 
cause for the apparent rotation of the nova shell between 
different epochs. Since this effect is generally not an issue 
with observations involving strong sources on larger scales 



< 
( 




Figure 6. Simulated radio images for the spherically symmetric Hubble flow model based on a 12 hour track. The contour levels are (-3, 
3, 3V2, 6, 6V2, 12, 12V2...) X 0.15 mjy. 




Figure 7. Simulated radio maps for a 24 hour track of the spherical Hubble flow model. The contour levels are (-3, 3, 3\/2, 6, &\/2, 12, 
1272...) X 0.15 mjy. 
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Figure 8. Simulated radio maps for the elliptical shell model based on a 12 hour track. The contour levels are (-3, 3, 3V2, 6, 6V2, 12, 
1272...) X 0.15 mjy. 
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Figure 9. Simulated radio maps for the elliptical shell model based on a 24 hour track. The contour levels are (-3, 3, 3V2, 6, 6V2, 12, 
1272...) X 0.15 mjy. 



12 Heywood & O'Brien 




Figure 10. Plots of the simulated uv coverage for a 12 hour 
(upper panel) and 24 hour (lower panel) track. 



Figure 11. Plots of the beam structure for the simulated 12 
hour (upper panel) and 24 hour (lower panel) observations. The 
beam structure is shown on the same scale as the simulated radio 
images. 



(e.g. radio galaxies) it is likely that these effects will only 
be present when observing objects with low surface bright- 
ness objects and small angular sizes, and would probably 
pass unnoticed when observing systems with less favourable 
dynamic timescales. Additionally a 24 hour track provides 
higher fidelity images due to more thorough sampling of the 
uv plane and an improved signal to noise ratio. 



4.2 Rings and discs 

The Fourier transform process can give rise to unexpected 
artefacts in the final image. In this section we briefly inves- 
tigate the possibility that the transition from a disc shaped 
structure to something that is ring-like could be responsible 
for unusual effects in the final radio images. Examination of 
Figure 3 shows how the model emitting region develops with 
time. Initially the source is optically thick and resembles a 
disc on the sky, with the flux density at this stage governed 
only by the angular size and the brightness temperature. The 
highest contour level also drifts around the field as the shell 
expands. Later a transition to optically thin emission occurs 
and the ring-like inner shell boundary is revealed. In order 



to assess the outcome of this transition on the final radio im- 
ages, simulated MERLIN observations of an expanding ring 
are presented for comparison with those of an expanding 
disc. These are presented in Figures 12 and 13. 

The model emitting region in each case has angular size 
scales and flux densities comparable with that of the shell 
of V723 Cas. As the intention here is to solely investigate 
the effect of the deconvolution process the final simulated 
data do not have Gaussian noise added. Note the differing 
contour levels for the ring and disc maps. The disc simula- 
tion is plotted with levels of (-3, 3, 3\/2, 6, 6\/2, 12, 12^/2...) 
x 0.15 mjy as with the real MERLIN maps and the sim- 
ulations presented in this paper. Since the simulated ring 
structure is much fainter the contour levels are reduced to 
(-3, 3, 3\/2, 6, 6\/2, 12, 12\/2...) x 0.03 mjy. The disc sim- 
ulation is brighter simply because there is more emitting 
material. 

Examination of Figures 12 and 13 reveals that even 
these simple zero noise models exhibit features such as ap- 
parent rotation and peaks dividing as the emitting region 
expands. These effects are seen in both the ring and the 
disc model meaning that it is unlikely that they appear in 
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Figure 12. Synthetic 12 hour MERLIN observations of an ex- 
panding ring. 
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Figure 14. Simulated uv coverage of a 12 hour track using the 
e-MERLIN 5 GHz contiguous band. 
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Figure 13. Synthetic 12 hour MERLIN observations of an ex- 
panding disc. 



the real observations due to the nova undergoing the transi- 
tion from disc-like to ring-like. This is further evidence that 
these effects are due to the limited uv coverage and are not 
related to the changing optical depth or other underlying 
properties of the source. 



5 CONCLUSIONS 

These results have shown that extreme care must be taken 
when using radio images alone to interpret the morphol- 
ogy of classical nova shells such as that of V723 Cas. The 
(comparatively) low radio surface brightness of a nova shell 



coupled with its low angular size gives makes it prone to in- 
strumental effects which may be mistaken for unusual shell 
morphology and motion (e.g. clumping, apparent rotation) . 
However these effects can be mitigated by using observations 
with longer tracks which both improves the uv coverage and 
increases the signal to noise ratio. When using radio imaging 
to track structure in novae and similar objects, coupling of 
simulations such as those presented here with observations 
would provide a more thorough insight into the underlying 
source structure. 

However, optical imaging of novae reveal that the emis- 
sion is far from homogeneous and exhibits equatorial and 
polar enhancements (e.g. FH Ser, Gill & O'Brien, 2000), 
clumping of material and even filamentary structures (e.g. 
GK Per, Slavin et al, 1995). Thus there is need for spherical 
models to be refined and it will be interesting to see how 
novae are observed by the next generation of interferome- 
ters such as e-MERLIN. With the vastly improved uv plane 
sampling of e-MERLIN exploiting multifrequency synthe- 
sis, Figure 14, one would not expect the instrumental effects 
documented in this paper to be present. Preliminary simula- 
tions shown in Figure 15 show three epochs of the elliptical 
shell model described in Section 2.3. Immediately apparent 
is the north-south aligned, extended, low brightness emis- 
sion which was below detectability levels in the simulations 
presented in Section 3. The root mean square background 
noise levels in the e-MERLIN simulations are ~5 fiJy. The 
uv plane sampling is consistent with that of e-MERLIN ob- 
serving using a 5 GHz contiguous band, i.e. 1000 frequency 
channels between 2 GHz and 5 GHz. The drawback here is 
that these simulations treat the source as being spectrally 
flat, which is obviously a poor assumption over such a wide 
frequency range. Nonetheless they show that the e-MERLIN 
view of sources of this nature will be much more faithful to 
the true source distribution due to its vastly improved uv 
coverage and sensitivity. 
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Figure 15. Preliminary e-MERLIN simulations of the December 
1996, December 1998 and October 2001 epochs of the elliptical 
shell model. 
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